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Q ; Abstract 



(N 



CN . An extension of the product operator formahsni of NMR is introduced, which 

uses the Hadamard matrix product to describe many simple spin 1/2 re- 
^ \ laxation processes. The utihty of this formahsm is iUustrated by deriving 

^^ ' NMR gradient-diffusion experiments to simulate several decoherence models 

C^ • of interest in quantum information processing, along with their Lindblad and 

On , Kraus representations. 

o 

O ' 1. Introduction 

;>: 

p ■ The product operator formalism is widely used for designing radio-frequency pulse se- 

quences to control the coherent evolution of multi-spin- 1/2 systems in liquid-state nuclear 
H ! magnetic resonance (NMR) spectroscopy ^^. This is a symbolic representation of the 

^1 2^ X 2^ Hermitian matrices of A^-spin density operators, observables and Hamiltonians, 

together with a collection of rules for evolving them under radio-frequency (RF) pulses and 
free evolution delays (in the Schrodinger, Heisenberg, or interaction picture defined by a 
^ ■ rotating frame Q]). This representation is obtained by expanding the matrices versus the 

"product operator" basis consisting of all possible A^-fold Kronecker (or tensor) products of 
the usual Pauli matrices cr^, cry, cr^ and the 2x2 identity cri. Together with the simpler and 
more general rules provided by the underlying geometric algebra structure p,|[ , the product 
operator formalism has also proven invaluable in designing pulse sequences to implement a 
wide variety of logic gates for quantum information processing (QIP) by NMR [[^-p!0[] (see 
rT| , p!2[| for recent reviews). 



O^. 
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In contrast, the general theory of NMR relaxation P,pr3| is usually developed using the 
"superoperator" representation of linear mappings on A^-spin density matrices p relative to 
the matrix basis \m){'n]/\ obtained from the eigenvectors of the Zeeman Hamiltonian \m) 
{■m,m/ G {0,1}-^). This representation consists of the 2^^ x 2^^ matrices acting on the 
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"coluinnized" density matrices versus the corresponding vector basis \'n][)\rri). The Lindblad 
differential operator |14] and Kraus operator sum |jl5|] representations, which are becoming 

m 



commonplace in quantum optics 



1^, and in foundational 



condensed matter physics 
studies of decoherence ||T8|-pO|, are virtually unknown to NMR spectroscopists. The recent 
use of NMR gradient-diffusion methods to simulate theoretical decoherence processes in ex- 
perimental studies of the control of decoherence for QIP ||2l| , |2^ has made it desirable to 
develop a new representation which: (i) blends naturally with the product operator formal- 
ism; (ii) avoids the unwieldy use of full 2^ x 2^ superoperator matrices; (iii) is sufficient 
to analyze common NMR gradient-diffusion experiments; (iv) describes the theoretical de- 
coherence models most often studied; (v) can easily be translated into the corresponding 
Lindblad, Kraus and superoperator representations. 

This paper presents such a formalism, which is based on the well-known Hadamard prod- 

This formalism permits the representation of relaxation processes 
superoperator matrix is diagonal relative to some operator {not 
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uct of matrices P3| , |24 
for which the 2^^ x 2^ 
superoperator) basis. As such it can handle both Ti and T2 relaxation (decoherence), but 
not cross-relaxation 0,|1^. In terms of Kraus operator sums, this is equivalent to requir- 
ing that the Kraus operators can all be simultaneously diagonalized. The utility of the 
Hadamard product formalism will be demonstrated by giving streamlined derivations of 
closed-form, time-dependent density operators for a variety of NMR gradient-diffusion ex- 
periments. These include collective and independent decoherence about any axis, together 
with several more complex forms of decoherence including the collective isotropic model. 
When possible, the corresponding Lindblad and Kraus operators for these experiments will 
also be given. We have found, however, that the Hadamard product permits a concise an- 
alytic description of decoherence with arbitrary correlations among the fluctuating fields at 
the different spins involved, whereas the corresponding Lindblad and Kraus forms entail the 
solution of polynomial equations, and can be written down explicitly only in the simplest 
cases. 



2. Hadamard products of product operators 

Given two M x M' complex-valued matrices A 



■iM-l.M'-l 
n*m'Jm,m'=0 

their Hadamard product is the matrix consisting of the complex products 

M-1,M'-1 



amm'lZ ^iZn aud B 

of corresponding pairs of elements, A B = [amm' &mm']mm'=o ~'^- ^^^^ the usual matrix 
product, the Hadamard product satisfies the mixed product identity with the Kronecker 
product ^,^ 



(A®B)0(C®D) 



(A0C)0(B0D) 



(1) 



but unlike the usual matrix product it is commutative. 

In the following, all matrices will be expressed as linear combinations of product oper- 
ators, as described above. Some further notations which prove useful include: M = 2^; 



CTl 



(T^ in the n-th place 



ai (the M X M identity); a 
n = 1, . 



CTi ■ • ■ (K 

N)- E^ = (I + <0/2 
E^ = (I -0/2 = 0-1 ••■ (|1)(1|) ■■■ 0cri (so 
EgE" = E"Eo = O, the M X M zero matrix); and 



^^,^■■ 


■ CTi (/i e {x, y, z} with 


= CTi ■ 


■■ (|0)(0) ■■■ 0cri, 


Kf = 


E^, (E^)2 = E^, and 



E, 



.U^m 



N 



TTll TT\N m 



5^ 



(2) 



where S^^ is the n-th bit in the binary representation m of the integer m, and the matrix 

<. On further abbreviating Eq by E = |00- • ■0)(00 ■ ■ -Ol, a 



^n\0 



powers [cr^ 
matrix A = [a. 



I, KY = 

mm'\m,m'=0 



iM-l 



]„ ^/=o ^^y ^^ expressed using these operations as 



A 



M-l 



m,m'=0 



Omm' CT— HjCT— 



(3) 



and in particular, a diagonal matrix D = Diag(d) (d = [d. 



lM-l> 
mm\m=0 j 



may be written as 



D 



E 



M-\ 



d E 



.A/-1 



mm -■-'m 



with d' = Hd, where H 

{not product) of all the spins p5 






•'— 'm=0 

[cr" + (T^')/\^) is the Hadamard transform 



The Hadamard product of any two such product operator expressions can be worked out 
from the mixed product formula (|l|) and the multiplication table below (in which ctq is the 
2x2 zero matrix). 

Table 1. 

Hadamard multiplication table 
for identity and Pauli matrices. 



© 


CTl 


CTx 


(Ty 


CTz 


CTl 


(Tl 


O-Q 


CTO 


CTz 


CTx 


o-o 


<Tx 


O-y 


CTO 


(Ty 


o-o 


O-y 


-o-^ 


<^0 


CTz 


O-z 


CTO 


f^o 


CTl 



The essential property of Hadamard products to be used in this Letter will now be given. Let 
A, B and C be M x M complex- valued matrices with A, C diagonal, and let a = diag(A), 
c = diag(C) be the column vectors formed from their diagonal elements. Then if c''" and C^ 
are the Hermitian conjugates of c and C, respectively: 



ABC^ = (ac^) ©B 



(5) 



If the matrices are sums of product operators as above, the dyadic product a c''^ may also be 
expressed in product operator form as 



ac' 



AHEHC^ . 



(6) 



3. Lindblad operators for gradient-diffusion 

Pulsed magnetic field gradients have many uses in NMR, for example to estimate con- 
ditional displacement probabilities in diverse transport phenomena p6[. The underlying 
theory is also important in relating NMR relaxation rates to the internal Brownian dy- 
namics of molecules |^. The experiments relevant to this Letter take advantage of the 



spatial extent of the ensemble of spin systems (molecules) in a liquid NMR sample. A mag- 
netic field gradient V-Bz parallel to the static field B^ along the z-axis causes the Zeeman 
precession rate of the spins to vary linearly with their spatial z-coordinates, winding the 
transverse (xy) magnetization into a spiral about z for which the net transverse magne- 
tization vanishes. The phase coherence thus rendered unobservable can be refocussed by 
either a second gradient pulse of the opposite polarity, or by an RF vr-pulse followed by a 
gradient of the same polarity. Hence to obtain true irreversible decoherence, it is necessary 
to wait for diffusion to randomize the positions of the molecules so that the correlation be- 
tween their spins' phases and z-coordinates is lost. A more detailed account of this process 
is outside the scope of this Letter, and may be found in Ref. p8 |. 



Let p be the M x M density matrix of an ensemble of spin systems each consisting 
of A^ spin 1/2 particles, assumed for simplicity to be noninteracting and to have the same 
gyromagnetic ratio 7, which has been polarized by a static magnetic field -B^ along the z- 
axis [^]. A uniform magnetic field gradient V-Bz correlates the phases of the spins with their 
spatial z-coordinates via the semiclassical propagator 

GW = y''~\mm{z)B^ = e-^^^^-l^-^-^)l\ (7) 



In this expression, ^ = 7 /J^ dr' dB^ j dz is the wave number of the phase ramp along the 



z 



N 



axis, gmm{z) = exp{—Lzk{N — 2h{m))/2) where h{m) = X]n=i ^m is the Hamming weight 
of m, and l"^ = —1. Since G{z) is diagonal, its action on a density operator p can be written 
as the Hadamard product G{z) pG'^i^z) = {g{z) g^{z)) p, and it is easily seen that 



g{z)g^z) = [e^^^^^^-^-'^^"^^J„,„,=o- (8) 

The effect of the molecular diffusion period t on the elements of g{z)g'^{z) is to convolute 
them with a Gaussian in z whose variance is Dt, where D is the diffusion constant p2| , p8| . 
These may be evaluated via Fourier transform methods to yield the corresponding phase 
damping matrix 

D(t) = [e-(^('^(-)-'^(-')»'^*]l';!=o • (9) 

Following a refocusing gradient pulse of equal magnitude and opposite polarity, the loss 
of coherence due to diffusion is now given (assuming no coherent evolution) by the time- 
dependent density matrix p(t) = D(t) p (so p(0) = p). Differentiation yields p{t) = 
— R p(t), where R = ['^mm']m,m'=o ^^ the rate matrix 

R = [(MM^)-M^'))f^]!!.io- (10) 



Conversely, integration yields D(t) = expQ(— Rt) = [exp(— rmm'^)]mm'=i • As t ^ oo all 
elements of D(t) and hence p{t) vanish save those with h{m) = h{m'). Those with m ^ m' 
represent coherences between states with equal angular momentum about the z-axis, which 
are called zero-quantum coherences in NMR [^. 

By expanding the square in Eq. (0), it is easily seen that R can be written as 

R = liiiieef + (iQe)!^) - ee'' , (ii) 

71 T -I 

where 1 is a length M vector of I's, i = [k y/2D h{m)]^^Q , and 1^, £ are their transposes. 
The application of Eq. (^ thus yields 

Pit) = Lp(t)L-lLVW-|p(i)L^ (12) 

which is the desired Lindblad master equation with a single real diagonal Lindblad operator 
L = Diag(£) = Lt [0,0,^. Because diag(cr^i + . . . + erf ) = [A^ - 2h{m)]^~^, this can 



also be expressed as L = k\J D 12 {N\ — (t\ — ■ ■ ■ — cr^), and because the action of L on 
p{t) is unchanged by its overall sign or by adding on a multiple of the identity I, it can be 
further simplified to 

L = k^Dj2{(Tl + ■ ■ ■ + O . (13) 

This result is readily generalized to cases in which the transverse magnetization from 
each type of spin has its own wave number fc". Then the propagator in (|^) becomes G{z) = 
exp{—Lz{k^(Tl + ■ ■ ■ + k'^cr^)), and an essentially identical derivation leads to Eq. (|1^) with 



L = ^/D/2{k'(T^ + ■■■ + fcV^) . (14) 

Decoherence processes of this kind occur naturally in heteronuclear gradient- diffusion exper- 
iments, and can be obtained in homonuclear by a sequence of gradient pulses interspersed 
with TT-pulses. If these refocus certain spins so that fc" = 0, the sequence applies the deco- 
herence selectively to the remaining spins. This analysis can be further generalized to the 
conditional gradient operations introduced in [^ , where the gradient pulses are interspersed 
with more complex RF pulse sequences and delays which implement conditional quantum 
logic gates such as the controUed-NOT |]5|. 

To illustrate such "conditional decoherence", consider a two-spin system subjected to 
a pair of gradient pulses selective for the first spin and interspersed by controUed-NOT's 
S^l^ = cr^^l + Eq to spin 1 conditional on 2, giving a net propagator [^ 

QU) = g-t2fe2<T,V2 gl|2g-t2:fci<T,V2 gl|2 

(15) 

where the subscripts on the /c's now specify the temporal order of the corresponding gradient 
pulses. A similar derivation then gives the Lindblad operator 

L = ^/DJ2(T^ik^I + k2(T^) (16) 



which for ki = k2 selectively decohers all off-diagonal elements of p except (t^EI ^^ pis = 
pIi- As another example, take a three-spin system and substitute the controlled-NOT's in 
Eq. (|15|) by Toffoli gates T^l^^ = I - E^Ef + cr^^EfE?. Then the Lindblad operator 



L = VD/2cT^ {kil + A;2(I + ct^ + cr ^ - cr^cr^)) (17) 

decohers the off-diagonal elements cr^EfEl ^^ p^j = pyg only if ki ^ k2, and otherwise 
creates this pseudo-pure state directly from p = a^. In general, the Lindblad operator 
will be proportional to the effective Hamiltonian of the gradient propagator preceding the 
diffusion period. 

Thus far the discussion has been restricted to a single diffusion period, so that the 
decoherence, although possibly selective or conditional, acts collectively on all the affected 
spins. The use of multiple diffusion periods permits implementation of arbitrary correlations 
including the independent case [^. The phase damping matrix in this latter case is the 
sum of those for each spin individually, leading to the Lindblad master equation 



.N 



Pit) = V (L„ Pit) L„ - iL„2 Pit) - Ipit) L„2) (18) 

•'— 'n=l 

with L„ = a/^/2 /c'^cr" for n = 1, . . . , A^. A general formula for arbitrary correlations is 
likewise expected to involve A^ Lindblads, but appears quite complicated. The result for a 
two-spin system, however, can be given as 

Li = aj^al + h^a^ , L2 = a^al + h+al , (19) 

with a±2 = li?i ± 5A, 6±2 = 1/^2 j^s/\, and 



R^R' - iSf 
4(5)' + (i?i - i?2)2 



(20) 



Here, EJ^ = (/c") D/2 are the selective decoherence rates for the two spins, and S with 
—VR^R^ < S < \'R}R? is a measure of the correlation in their mutual decoherence. These 
equations are indeterminate ii R^ — R^ = Q = S, and care must be taken in the choice of 
signs for the square roots of a^^ and h±^. Specifically, if R^ > R^ and S'^ > R^iR^ — R^)/2 
all four roots are positive, whereas 6_ is negative if S"^ < R^iR^ — R^)/2, and if R^ < R^ 
then all four roots are positive unless 5*^ < R^iR'^ — R^)/2, in which case the negative square 
root is taken for a_ ; if 5* = then Eq. (|TH|) is used instead. 

In closing this section, we note that standard results in the theory of Kronecker products 
|2^ show that the M^ x M^ matrix C oi a. Lindbladian superoperator of the form given in 



Eq. (|T8|) (with any Hermitian L„) is 



-r X \^ / -T- \ 9 

^ = iE jLl®I-I®Ln , (21) 



and that this is diagonal whenever the L„ are. 



4. Kraus operator sums and correlated decoherence 

Performing an eigenvalue decomposition of the phase damping matrix (^) gives the in- 
tegrated evolution equation in the standard Kraus operator sum form, 

pit) = pQ D(t) = p V ^ k„(t) ^m{t) klit) 

(22) 

yM-l 

K„(t)pK„(t) , 

'm=0 






where the real diagonal matrices Km(t) = ^JKm{t) Diag(km(t)) are in general complicated 
functions of time t. The only assumption made here is that D(t) is positive-semidefinite for 
all t, as it must if p{t) is to be positive-semidefinite for all t and inital states p = p(0). 

Due to its algebraic complexity, in general the Kraus form can be obtained only by 
numerically diagonalizing D(t) at each time point. For independent decoherence at the 
different spins, however, the rates rmm' are easily seen to be proportional to the squared 
Hamming distances J2n=ii^m ~ ^m')^5 rather than to {h{m) — h{m')Y. This leads to the 



Kraus form 29 



.jv^-*^-i 



p{t) = 2~^}_^ ^Ut) <P < , (23) 

where Kmit) = Y\n=i{^ + (— l)'^™p"(i)) for the one-spin survival probabilities p"(t) = 
exp(— ti?"). The case of collective decoherence appears substantially more difficult, and 
the eigenvalues i^m{t) involve radicals even for just two spins. Assuming as in Eq. (P) that 
the one-spin survival probabilities are both equal io p = p{t), a simple algebraic form, not 
based on diagonalization, exists and is given by 

Pit) = liil+p)I- pa^<T^)piil+p)I- pa^a^) 

+ 1(1 - ip)") (I + ^>^) P (I + 'T^z') (24) 

This does not seem to extend to larger numbers of spins. 
The extended Kraus (diagonal) operator sum form, 

P(^) = E , c^^,it)ai^pa^ , (25) 

turns out to be algebraically simpler, and (given that [cmm'(^)] is positive-semidefinite) is 
readily converted into the above standard form by diagonalization. To derive the extended 
form with arbitrary correlations between two spins, Eq. (^ is used to express the microscopic 
effect of an arbitrary decoherence process as 



4 






(26) 



where the exponentials on the right-hand side have been commuted past each term of 
HEH = il + cr^ + cr^ + (t^(t^)/4:, changing the sign of their arguments as appropriate. 
On combining the pulse sequence which implements this propagator with one or more dif- 
fusion periods, the resulting macroscopic average is given by ||22[ 



pit) = i (l + e-^'* al + e-«'* a^ + e-(«^+«^+^-^'-^')* o-.V^j p , (27) 

where the rates R^, R^ and S are as given previously. 

The Hadamard multiphcation table (Table 1) shows the identity term in ( p7|) annihilates 
all but the I, cr^, cr^ and cr^cr^ components of p, while the cr^ term annihilates all but the 
cr^, (T^, cr^cr^ and cr^cr^ components, etc. These components are isolated by the projections 
parallel (+) and perpendicular (— ) to the z-axis, 

P,,e2 = KP + ei'^zVcrz' + e2cr>cr2 + eie2CT-»cr>2) (28) 

(ei,e2 G {±1}), which enables the evolution equation to be written without using the 
Hadamard product as 

p(t) = P++ + e-^'* p^ + e~«'* p+_ + e-^'* e-«'* e"''''''"' p . (29) 

Remarkably, Eqs. (p7|^9|) extend easily to any number of spins. The last term in 



(p7|), previously derived in p2[, can be cast in a more symmetric form by expanding 

exp{—S(T^(T^t) = cosh(5't) — cr^^cr^^sinh(S't) and applying it to p in Eq. (^), leading 

eventually to the two-spin extended Kraus operator sum 



8p{t) = a {(tIpo-I + (T^p ctI - ctIctIp - p a^a^) + 6++ p 



(30) 



where a = pV(l - {qf)/Q and b^^^ = 2(1 + eip^)(l + eap^) + e^p^p^iq - if/q (ei, e2, eg G 
{±1}), with p^ = exp(-i?H), p2 = exp{-R^t), q = exp{-St). 

Finally, by sandwiching a gradient between RF pulses for a rotation and its inverse, the 
spins can be decohered about any axis. For example, a Hadamard transform may be used 
to obtain collective decoherence about the x-axis as 

Pit) = H(D(t)0(HpH))H, (31) 

where D(t) is as given in Eq. (|). Isotropic collective decoherence about all three axes is 
obtained by applying three identical gradients, two of which are sandwiched between RF 
pulses for (7r/2)-rotations about x and y, and each followed by equal diffusion periods. The 
Lindblad operators in this case are that in Eq. ([13D together with the two obtained by 
replacing cr^ by (Ty^ and cTy throughout. The integrated form may be written compactly as 

pit) = ZHZ^(D(t)0(ZHZ"fH(D(t)0(H(D(t)0p)H)HZHZ"f)))ZHZ^ (32) 
where Z = exp(-i7r(cr^^ H h (t^)/A)). 

5. Conclusions 

It has been shown that NMR gradient-diffusion methods enable precise implementa- 
tions of the adiabatic decoherence processes most often studied in QIP [^, and that the 
Hadamard product formalism is a simple and efficient means of analyzing such processes 
(regardless of their underlying physical mechanism) . Nonadiabatic relaxation of p towards 

8 



the equilibrium density matrix p^^ can also be described using the Hadamard product, 
simply by decohering the diagonal part of p — p about the x-axis. Using the fact that 
I X = Diag(diag(X)) for any X, this leads to 

Pit) = H(D(t) (H(I (p - p,q))H))H + p,q , (33) 



rN 



where D(t) = fl^^-^ expQ(— t(I + cr")/T") for independent decoherence (although any of the 
forgoing decoherence models could be used). Nonadiabatic T2 processes can similarly be 
described using complex- valued phase damping matrices D(t). 

The ability to simulate complicated decoherence processes by gradient-diffusion, together 
with modulation of the natural decoherence processes operative in NMR |^, should sig- 



nificantly enhance the utility of NMR as a testbed for the development of more power- 



ful quantum information processors |12[. As a means of designing such experiments, the 
Hadamard formalism is effectively limited to operations for which the propagator G{z) can 
be diagonalized by a unitary matrix that is independent of z, since otherwise the algebra 
and integrations will not usually be analytically tractable. Some decoherence processes, 
e.g. those involving cross-relaxation P,p!3|, cannot be efficiently described by Hadamard 
products, and it is not known if they can be obtained by gradient-diffusion. Further work 
on the universality [^^, and complexity, of these simulations is needed. 
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